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Abstract 

We give a survey of recent results on flow-structure interactions modeled by a modified 
wave equation coupled at an interface with equations of nonlinear elasticity. Both subsonic 
and supersonic flow velocities are considered. The focus of the discussion here is on the 
interesting mathematical aspects of physical phenomena occurring in aeroelasticity, such 
as flutter and divergence. This leads to a partial differential equation (PDE) treatment of 
issues such as well-posedness of finite energy solutions, and long-time (asymptotic) behavior. 
The latter includes theory of asymptotic stability, convergence to equilibria, and to global 
attracting sets. We complete the discussion with several well known observations and 
conjectures based on experimental/numerical studies. 

Key terms: flow-structure interaction, nonlinear plate, flutter, well-posedness, long-time 
dynamics, global attractors. 


1 Introduction 

The main goal of this paper is to provide a mathematical survey of various PDE models for 
fluttering plates The flutter pheonmenon is of great interest across many 

fields. Flutter is a sustained, systemic instability which occurs as a feedback coupling between a 
thin structure and an inviscid fluid. A static bifurcation may also occur, known as divergence or 
buckling. Here we consider flow-plate dynamics corresponding to large range of flow velocities 
including both subsonic and supersonic. 

Flow-structure models have attracted considerable attention in the past mathematical lit¬ 
erature, see, e.g., and the references therein. However, the 

majority of the work that has been done on flow-structure interactions has been devoted to 
numerical and experimental studies, see, for instance, 

HHm El S3 and also the survey [25] 
and the literature cited there. Many mathematical studies have been based on linear, two di¬ 
mensional, plate models with specific geometries, where the primary goal was to determine the 
speed at which flutter occurs see also 0IH 0 0 El for the recent studies of 

linear models with a one dimensional flag-type structure. 

Given the difficulty of modeling coupled PDEs at an interface [50l|5l, theoretical results 
have been sparse. While numerical studies address various situations, they are based on finite 
dimensional approximations of continuum models fully described by PDEs 0 0 HDJ. The 
infinite dimensional nature of the physical phenomena may not be adequately reflected by these 
approximations. The results presented herein have demonstrated that the models can be studied 
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from an infinite dimensional point of view, and moreover that meaningful statements can be 
made about the physical mechanisms in flow-structure interactions strictly from the PDE model. 
The problems in the PDE analysis revolve around (i) the mismatch of regularity between two 
types of dynamics: the flow and the structure which are coupled in a hybrid way, (ii) the 
physically required presence of unbounded or ill-defined terms (boundary traces in the coupling 
conditions), and (iii) intrinsically non-dissipative generators of dynamics, even in the linear case. 
(The latter are associated with potential chaotic behavior.) 

This paper specifically addresses the interactive dynamics between a nonlinear plate and a 
surrounding potential flow mm- The class of models under consideration is standard in the 
modeling of flow-structure interactions, and goes back to classical discussions (13 [34], and also 
[35) 41 and the references therein). The description of physical phenomena such as flutter, 
divergence or buckling will translate into mathematical questions related to existence of nonlin¬ 
ear semigroups representing a given dynamical system, asymptotic stability of trajectories, and 
convergence to equilibria or to compact attracting sets. Interestingly enough, different model 
descriptions (including boundary conditions) lead to an array of diverse mathematical issues 
that involve not only classical PDEs, but subtle questions in non-smooth elliptic theory, har¬ 
monic analysis, and singular operator theory. Thus, in this paper, we concentrate primarily on 
mathematically rigorous aspects of modeling aeroelastic dynamics and their analysis. For more 
details concerning the relationship of the theory developed with numerics and experiments we 
refer to the paper [Si- 

Notation: For the remainder of the text denote x for (x,y,z) £ or ( x,y ) £ fi C 
as dictated by context. Norms || • || are taken to be L 2 (D) for the domain D. The symbol v 
will be used to denote the unit normal vector to a given domain, again, dictated by context. 
Inner products in L 2 (R+) are written (-,-), while inner products in ^(R 2 = <9R(j_) are written 
(-,-). Also, H S (D) will denote the Sobolev space of order s, defined on a domain D , and 
Hq(D) denotes the closure of Cq°(D) in the H S (D) norm denoted by || • ||//«(d) or || • || Si d. We 
make use of the standard notation for the boundary trace of functions, e.g., for (f> £ ff 1 (R')_), 
7 (</>) = tr[4>] = 4>\, =0 is the trace of on the plane {x : 2 = 0}. 

2 PDE Description of the Gas-Plate Model 

To describe the behavior of the gas (inviscid fluid), the theory of potential flows [9) (34j 01] is 
utilized. The dynamics of the plate are governed by plate equations with the von Karman (vK) 
nonlinearity EH EH- This model is appropriate for plate dynamics with “large” displacements, 
and therefore applicable for the flexible structures of interest [9j EH El 0TJ09]. 

The environment considered is = {(x,y,z) £ R 3 : z > 0}. The thin plate has negligible 
thickness in the z-direction. The unperturbed flow held moves in the ^-direction at the fixed 
velocity U. The physical constants have been normalized (corresponding to the linear potential 
flow) such that U = 1 corresponds to Mach 1. The plate is modeled by a bounded domain 
C <9R(j_ with smooth boundary <9f l = T. 

The scalar function u : fix R + —> R represents the displacement of the plate in the z-direction 
at the point (x, y) at the moment t. The plate is taken with general boundary conditions BC{u) 
which will be specified later. The dynamics are thus given by: 

f (1 - aA)uu + A 2 x y u + D 0 (u t ) + f(u) = p(x, t) in D 

\BC(u) on dn, ' 

with appropriate initial data uq, u\. The quantity p(x, t) corresponds to the aerodynamic 
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pressure of the flow on the plate and, in this configuration, is given in terms of the flow: 

p(x,t) =p 0 (x) + (d t + Ud x )tr[(j)\, x £ fl. (2.2) 

The quantity p o represents static pressure applied on the surface of the plate. The function 
D 0 represents an abstract form of interior damping. A parameter a > 0 represents rotational 
inertia in the filaments of the plate and is proportional to the square of the thickness of the 
plate. The case of interest, a = 0, is said to be non-rotational ; the mathematical effects of a 
are discussed below. The vK nonlinearity is given by: 


f(u) = ~[u,v(u) +-F 0 ], 

where F 0 is a given in-plane load of sufficient regularity. The vK bracket is given by 

[u, w) 'a xX Wyy “1“ UyyW XX 2U X yW X y, 

and the Airy stress function u(ui,u 2 ) solves the elliptic problem 

A 1 2 u(ui,u 2 ) + [u\,U 2 \ = 0 in f l, d u v(u\ 7 u 2 ) = v(ui, u 2 ) = 0 on T, 

(the notation v(u) = v(u,u) is employed). Here and in In what follows v (resp., r) denotes the 
unit outer normal (resp., tangential) to H. 

The following types of boundary conditions for the displacement u are typical for flow- 
structure interaction models. 

(C) Clamped boundary conditions (corresponding to a panel element in the flow-structure in¬ 
teraction) take the form 

u = d v u = 0 on T. 

In this case, the plate is considered to be embedded in or affixed to a large rigid body. 

(FC) Let the boundary be partitioned in two pieces: To and Ti. Free-clamped boundary condi¬ 
tions (possibly corresponding to a flap, flag, or cantilevered airfoil) are given by 

! u = d v = 0 on To 
Au + (1 — p)Biu = D\(d v u t ) on Ti, 
d u Au + (1 — h)B 2 u - piu = D 2 (u 7 u t ) + ad u u t t on Ti, 

where the boundary operators B\ and _B 2 are given by [49] : 

Biu = 2v\V2U xy — ufu yy — v% u xx = —d x u — V • v(x)d v u , 

(2.3) 

B 2 U = d T [(^1 ^ 2 ) T ^1^2 (l^yy W^x)] — d x di/d x U] 

The parameter is nonnegative and 0 < p < 1 is the Poisson modulus. The operators 
Di for i = 1,2 are possible energy damping/dissipation functions. 

The scalar function 4> ■ R+ x M —>■ R is the linear flow potential and satisfies: 

f ( d t + Ud x ) 2 (j) = A Xt y tZ (j) in R 3 + . 

\FC(f) 7 [ ' J 

with appropriate initial data (f> 0 and ct>E The term FC((f>) represents the flow boundary condi¬ 
tions or interface conditions. We will consider two primary flow conditions: 

1 Note that, as mentioned above, U is normalized by the speed of sound and thus is what engineers refer to 

as the Mach number. 
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(NC) The standard flow boundary conditions—the full Neumann condition —are henceforth de¬ 
noted (NC). This is typically utilized when a majority of the plate boundary is clamped 
or hinged. It takes the form: 

dz<p\ z=0 = ( dt + Ud x )u ex t, on R 2 = <9R+, (2.5) 

where the subscript ‘ext’ indicates the extension by zero for functions defined on to all 
of R 2 . 

(KJC) The Kutta-Joukowsky flow condition, henceforth denoted (KJC), is a mixed (Zaremba-type 
EDI) condition which will be discussed at length below. We introduce the acceleration 
potential ip = (dt + Ud x )(p for the flow potential (p . Then the Kutta-Joukowsky flow 
condition on R 2 is 


d z (p = (d t + Ud x )u on II; ip = 0 on R 2 \ Q. (2.6) 

In both cases above the interface coupling on the surface of the plate occurs in a Neumann type 
boundary condition known as the downwash of the flow. The choice of flow conditions itself 
dependent upon the plate boundary conditions, and this is determined by application. 


3 Standard Panel Model 

Our primary interest here is the following PDE systerr@ 

'(1 - aA)u tt + A 2 u + ku t + f(u) = p 0 + (d t + Ud x )tr[<p] 
u(0) = u 0 ; Ut(0) = ui, 
u = d v u = 0 
(d t + Ud x ) 2 cp = A (p 
(p( 0) = </>o; <Pt( 0) = (pi 

d v (p = -[(<% + Ud x )u(x)\ ext 

In the discussion below, we will encounter strong (classical), generalized (mild), and weak (vari¬ 
ational) solutions. In the results below, semigroup theory is utilized, which requires the use 
of generalized solutions; these are strong limits of strong solutions. These solutions satisfy an 
integral formulation of (13.lj> . and are called mild by some authors. We note that generalized 
solutions are also weak solutions, see, e.g., ;22 s Section 6.5.5] and [73]. For the exact defini¬ 
tions of strong and generalized solutions we refer to mm Hi [73]. Finite energy solutions are 
identified with mild solutions. 

3.1 Well-Posedness of Nonlinear Panel Model 

Well-posedness results in the past literature deal mainly with the dynamics possessing some reg¬ 
ularizing effects. This has been accomplished by either accounting for non-negligible rotational 
inertia mmm and strong damping of the form —aAu t , or by incorporating helpful thermal 
effects into the structural model EH EH]- In the cases listed above, the natural structure of 
the dynamics dictate that the plate velocity has the property ut € if 1 (12), which provides the 
needed regularity for the applicability of many standard tools in nonlinear analysis. One is still 

2 For definiteness and some simplification we concentrate on the clamped boundary conditions for the dis¬ 
placement u. 


ill 12 X (0, T), 

on <912 x (0,T), 
in R)j_ x (0, T), 
in R)j_ 

011 ®{(x,y)} x (0) T)- 
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faced with the low regularity of boundary traces, due to the failure of the Lopatinski conditions 

[55] . In fact, the first contribution to the mathematical analysis of the problem is mm (see 
also [22[ Section 6.6]), where the case a > 0 (rotational) is fully treated. The method employed 
in mmm relies on sharp microlocal estimates for the solution to the wave equation driven 
by U 1 / 2 ( Cl) Neumann boundary data given by u t + Uu x . This gives <j> t |n G L 2 (0, T; _ff _1 ' /2 (n)) 

[56] . Along with an explicit solver for the three dimensional flow equation (the Kirchhoff for¬ 
mula), and a Galerkin approximation for the structure, one may construct a fixed point for the 
appropriate solution map. Such a method works well for all values of U. When a = 0, and 
additional smoothing is not present, the corresponding estimates become singular, destroying 
the applicability of the previous methods. The main mathematical difficulty of this problem 
is the presence of the boundary terms: (< fit + U<j> x ) |n acting as the aerodynamic pressure on 
the plate. When U = 0, the corresponding boundary terms exhibit monotone behavior with 
respect to the natural energy inner product (see [22j Section 6.2] and [5(J]) which is topologically 
equivalent to the topology of the state space Y. The latter enables the use of monotone operator 
theory [22] (and references therein). However, when U > 0 this is no longer true with respect 
to the topology induced by the energy spaces. 

We begin with an overview of the well-posedness results for the flow-plate model which covers 
both subsonic and supersonic cases. 

Theorem 3.1. Consider problem in \3.1\) with U / 1, p € L 2 (VL),Fq 6 H 3+S (fl). Let T > 0 
and 

wo, Ui) G Y = Yfi x Y pl = ( H\R 3 + ) x L 2 (R 3 + )) x (H 2 (H) x L“(H)) (3.2) 

where Lfikl) = Hq(LI) if a > 0 and = L 2 (Q) when a = 0. Then there exists unique 

generalized solution 

(cf(t),Mt)-,u(t),u t (t)) G C([0,T],Y). (3.3) 

This solution is also weak and generates a nonlinear continuous semigroup St : Y Y. 

For any initial data in 


Y = 


jy = (<Mi;u,ui) G Y 


men, 

-U 2 d x 2 ^ + A^eL 2 (R 3 + ), 
d u f> = -[mi + Ud x u] ■ In G H 1 (K 2 ), 
-A 2 u + Utr[d x (f\ e [L£(ft)]* 


(3.4) 


the corresponding solution is also strong. Here [L^n)]* is H 1 (f2) if a > 0 and L 2 (kl) when 
a = 0. 


In the subsonic case, there is no spatial “degeneracy” in the flow equation, leading to a 
stronger well-posedness result: 

Theorem 3.2. Suppose U < 1, p G L 2 (Q) and Fo G H 3+S (Q). 

Generalized and Weak Solutions: Assume that initial data sytisfy (13.21) . Then (EZED has 
a unique generalized solution on any interval [0,T]. Every generalized solution is also weak. 
Moreover these solutions to <ED generates a dynamical system ( St,Y ). Any generalized (and 
hence weak) solution to (13.11) . satisfies the bound 

sup {((1 - aA)u u u t ) o i0 + II Au||n + \\(ft\\li + II'W’lljU } < +oo. (3.5) 

t>o t + + J 

Strong Solutions: Assume u\ G Hq(LI), uq G W = {w G H 2 (H) : A 2 w G L 2 (fA)}. Also, 
denote H = Hq(LI) x LflfA). Moroever, suppose <f >o G H 2 (M^) and <j >i G H 1 (M'^), with the 
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following compatibility condition in place: 


<9 2 0o 


z=0 


u\ + Ud x uo if x £ fl 
0 if x ^ n 


(3.6) 


Then (13.11) has a unique strong solution for any interval [0,T], This solution possesses the 
properties 


(<j>\ <Pu4>tt) S Lco^OjT', H~ 


l 3 + ) x H 1 ) 




) x L 2 (K+)), (u;u t ;u tt ) £ Too(0,T;IT x U), 


and satisfies the energy identity 


£{t) + 



£(s), t> s, 


where Eft) = E p i(t) + Efi(t) + E int (t), with: 


(3.7) 


e pI = j[((l-aA)M fl «,)+ ||Au||q n + i||Aw(u)||^ n ] - (F 0 ,[u,u])n + {p,u)n (3.8) 
E fl = ^[II^II^IRS +||V0|| o , K 3 -C/ 110CC11 0 ,]R^ ]! E'int = 2U(<p, u x )ci- (3-9) 

The proof of Theorem 13.21 given in [22] ans m relies on three principal ingredients: (1) 
renormalization of the state space which yields shifted dissipativity for the dynamics operator 
(which is nondissipative in the standard norm on the state space); (2) the sharp regularity of 
Airy’s stress function, which converts a supercritical nonlinearity into a critical one [22] (and 
references therein): 

||?;(u)|| M /2,oo( f 2) < C\\u\\%i( n y, (3.10) 

and (3) control of low frequencies for the system by the nonlinear source, represented by the 
inequality 

IMl 2 (fi) ^ e [ll u llff 2 (n) 2 + l|Aw(w)||i 2 (fi)] + Ce- (3-11) 

An alternative proof based on a viscosity method can be found in [2 3) . 

In comparing the results obtained for supersonic and subsonic cases, there are two major 
differences at the qualitative level, in addition to an obvious fact that Efi is no longer positive 
when U > 1: First , the regularity of strong solutions obtained in the subsonic case [23||13] 
coincides with regularity expected for classical solutions. In the supersonic case, there is a 
loss of differentiability in the flow in the tangential x direction, which then propagates to the 
loss of differentiability in the structural variable u. As a consequence strong solutions do not 
exhibit sufficient regularity in order to perform the needed calculations. To cope with the 
problem, special regularization procedure was introduced in [26] , where strong solutions are 
approximated by sufficiently regular functions, though not solutions to the given PDE. The 
limit passage allows one to obtain the needed estimates valid for the original solutions [26]. The 
analysis of the supersonic case makes use of a change of state variable: the dynamic variable 
if = 4>t + Ucj) x is considered as the second flow state. In this case the adapted energies are: 


E fi — 2[II^Ho,R]_ + ^ ~ E pi~h E f£ 

and the energy identity is given by: 

£{t)+ [ [ ku 2 t + f {u x , ip)ndr = £{0). 

Js JQ JO 


(3.12) 


(3.13) 
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Secondly , in the subsonic case one shows that the solutions are bounded in time, see (22l 
Proposition 6.5.7] and also [231173] . This property cannot be shown in the supersonic analysis. 
The leak of the energy in energy relation can not be compensated for by the nonlinear terms 
(unlike in the subsonic case). 

3.2 Long-time Behavior and Attracting Sets 

It is known [41] (and references therein) that, although the flow coupling may introduce potential 
instability (flutter), it can also help bound the plate dynamics and may assist in the dissipation of 
plate energy associated to high frequencies J3S] [37] □ The work in [34] is an attempt to address 
these observations rigorously, via the reduction result in Theorem 13.61 below. This result is 
implemented and the problem is cast within the realm of PDE with delay (and hence our state 
space for the plate contains a delay component which encapsulates the flow contribution to 
the dynamics). One would like to show that the model described above is stable in the sense 
that trajectories converge asymptotically to a “nice” set. Unlike parabolic dynamics, there is 
no a priori reason to expect that hyperbolic-like dynamics can be asymptotically reduced to 
truly finite dimensional dynamics. By showing that the PDE dynamics converges to a finite 
dimensional (compact) attractor it effectively allows the reduction of the analysis of the infinite 
dimensional, hyperbolic-like, unstable model (asymptotically in time) to a finite dimensional set 
upon which classical control theory can be applied. 

Remark 3.1. Based on numerical results, it is well known that for a panel whose width is much 
larger than its chord, that four to six structural modes are sufficient to give very good accuracy 
to predict flutter and limit cycle oscillations. Indeed, for transonic and subsonic flows even one 
or two modes may be sufficient. As the width of the panel is decreased more and more modes 
are required however. Also if the plate is under tension, then the number of modes needed 
increases, and in the singular limit of infinite tension the number of modes required also goes 
to infinity [35] [40] . 

The main point of the treatment in [24] is to demonstrate the existence and finite-dimensiona¬ 
lity of a global attractor for the reduced plate dynamics (see Theorem 13.61) in the absence of 
any imposed damping mechanism and in the absence of any regularizing mechanism generated 
by the dynamics itself. Imposing no mechanical damping, taking a = 0, and assuming the flow 
data are compactly supported yields the primary theorem in [24] : 

Theorem 3.3. Suppose 0 < U 7 ^ 1, Fo £ H 3+S (kl) and po £ L 2 (fl). Then there exists a 
compact set C Y p i of finite fractal dimension such that 

lim d Ypl ({u{t),u t (t)),^) = lim inf (||n(f) - i^olli + IK(*) - ^i|| 2 ) = 0 

t— >00 * t -¥OO 

for any weak solution (u,ut; <j>, 4>t) to (ED with initial data 

Y = H 2 (n) x L 2 (fi) x h\r 3 + ) x l 2 (r 3 + ), 

which are localized in K+ (i.e., </>o(x) = cj> i(x) = 0 for |x| > R for some R > 0). Addition¬ 
ally,there is the extra regularity C (i7 4 (D) nHg(fl)) x H 2 (fl). 

The proof of Theorem 13.31 requires modern tools and new long-time behavior technologies 
applied within this delay framework. Specifically, the approach mentioned above (and utilized 
in [HED [58]) does not apply in this case. A relatively new technique Cm ED ED] allows one 

3 Though, there are also calculations performed at the physical level which demonstrate that the effect of the 

flow can be destabilizing in certain regimes—see Remark 14. II below. 
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to address the asymptotic compactness property for the associated dynamical system without 
making reference to any gradient structure of the dynamics (not available in this model, owing to 
the dispersive flow term). In addition, extra regularity and finite dimensionality of the attractor 
is demonstrated via a quasi-stability approach |22| . The criticality of the nonlinearity and the 
lack of gradient structure prevents one from using a powerful technique of backward smoothness 
of trajectories |22| . where smoothness is propagated backward from the equilibria. Without a 
gradient structure, the attractor may have complicated structure (not being characterized by 
the equilibria points). In order to cope with this issue, a novel method that is based on suitable 
construction of an e net and exploits only the compactness of the attractor. 

In the presence of additional damping imposed on the structure, it is reasonable to expect 
that the entire evolution (both plate and flow) is strongly stable in the sense of strong conver¬ 
gence to equilibrium states. Such results are shown for the model above (in the same references) 
only for the case of subsonic flows , U < 1. The key to obtaining such results is a viable energy 
relation and a priori bounds on solutions which yield finiteness of the dissipation integral. We 
obtain a result for stabilization of the full flow-plate dynamics (as opposed to just the structural 
dynamics, as above, for all flow velocities U ^ 1). For this we recall the global bound in (13.51) 
which along with the energy identity in (13.71) allows us to obtain the following: 

Corollary 3.4. Let k > 0 and a = 0. Then the dissipation integral is finite. Namely, for a 
generalized solution to en we have 

WMtfWo&dt < K u>k < oo. 

From this finiteness, we can show §53]: 

Theorem 3.5. Let 0 < U < 1 and a = 0. Assume p 0 £ L 2 (^) and F 0 £ Then for all 

k > 0, any solution {u(t), 4>(t), 4>t(t)) to the flow-plate system (13.11) with initial data 

{u 0 ,u i;0 o ,0i) £ (H$nH 4 ){n) x H 2 {Q) x H 2 (R 3 + ) x H\R 3 + ). 

that are spatially localized in the flow component (i.e., there exists a po > 0 so that for |x| > po 
we have </>o(x) = (fifx .) =0,) has the property that 

+ H m ‘WIIL(0) + H(t) ~ (k p ) + \\Mt)\\l 2 ( Kf ,)} = 0 

for any p > 0, where K p = {x £ : ||x|| < p] and M denotes the set of stationary solutions 

to (13.11) (for their existence and properies see 

The above result remains true for finite energy initial data (only in Y) with convergence in a 
weak sense. Additionally, for the system with rotational inertia present a > 0 (and correspond¬ 
ing strong damping) the analogous strong stabilization result holds QjJ [22] for finite energy 
initial data 

{u 0 , Ul ^oHi) € Hom X L£(fi) X WflRl) x L 2 (R 3 + ). 

The fact that k > 0 is needed in order to obtain convergence to equilibria is corroborated 
by the counterexample EE, which shows that with k = 0 periodic solutions may remain in the 
limiting dynamics. In fact, k > 0 allows one to prove that the entire flow-structure system is a 
gradient system. In view of the above result—Theorem 13.51 -with finite energy initial data, we 
see that any damping imposed on the structure seems to eliminate the flutter. This is consistent 
with experiment for clamped plates. From a physical point of view, the aerodynamic damping in 
subsonic flow for a plate is much lower than in supersonic flow, so the plate tends to oscillate in 



a neutrally stable state in the absence of aerodynamic damping or structural damping j40j. On 
the other hand, when U > 1 we expect that such result no longer holds (see, e.g., the numerical 
analysis in [5911. 

We also note that when the vK nonlinearity is replaced by Berger’s nonlinearity, and the 
damping coefficient k > 0 is taken sufficiently large, it has been shown |54| that all finite energy 
trajectories converge strongly to equilibria, and hence the flutter is eliminated. Though we 
have Theorem 13.51 for smooth initial data , the analogous result for the vK nonlinearity is still 
unknown, unless additional static damping controller is used [551 . 

3.3 Reduced Delayed Model 

A key to obtaining attracting sets is the representation of the flow on the structure via a delay 
potential (see Section 3.3 in [22]). Reducing this full flow-plate problem to a delayed von Karman 
plate is the primary motivation for our main result and permits a starting-point for long-time 
behavior analysis of the flow-plate system, which is considerably more difficult otherwise. The 
exact statement of this key reduction is given in the following assertion: 

Theorem 3.6. Let the hypotheses of Theorem \3.3\ be in force, and (uo, up, fo, </>i) £ H x L 2 (fi) x 
77 1 (K^_) x Zy 2 (®M-)- Assume that there exists an R such that <^o( x ) = <^i(x) = 0 for |x| > R. 
Then the there exists a time t# (R,U,Q) > 0 such that for all t > the weak solution u(t) to 
mu satisfies the following equation: 

u tt + A 2 u - [u , v{u) + T 0 ] =po- (d t + Ud x )u - q u (t) (3-14) 

with 

1 f^ /*2 7T 

q u (t) = — / ds dd[MgU ext ](x — (U + sm9)s,y — scosd, t — s). (3.15) 

2tt Jo Jo 

Here, Mg = sin 6d x + cos 9d y and 

t* = inf{f : x(/7, 9, s) Q for all x £ fi, 9 £ [0, 2tt], and s > t} (3.16) 

with x(t/, 9, s') — {x — {JJ T sin#)s, y — s cos 9) C K 2 . 

The flow state variables ( (f>,(f>t ) manifest themselves in our rewritten system via the delayed 
character of the problem; they appear in the initial data for the delayed component of the 
plate, namely u t \^_ t , (J ,. The reduced model displays the following features: (i) it does not 
have gradient structure (due to dispersive and delay terms), (ii) the delay term appears at the 
critical level of regularity. This raises a natural question regarding long time behavior and the 
existence of global attractors. However, despite of the lack of gradient structure, compensated 
compactness methods allow to “harvest” some compactness properties from the reduction, so 
the ultimate dynamics does admit global attracting set. In fact, the result of Theorem 13.31 is a 
consequence of a more general result pertinent to delay PDE system, and formulated in Theorem 

m 

4 Simplified/Reduced Models 

4.1 Plate Model with Delay as a Dynamical System 

Below we utilize a positive parameter 0 < t* < +oo as the time of delay, and accept the 
commonly used (see, e.g., [33] or m) notation u t {-) for function on s £ [—1*,0] of the form 
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s i ^ u(t+s). This is necessary due to the delayed character of the problem which requires initial 
data of the prehistory interval [—i*, 0], i.e., we need to impose an initial condition of the form 
u| tg (_ t . 0 ) = ? 7 (x,t), where rj is a given function on SI x [—t*,0]. We can choose this prehistory 
data rj in various classes. In our problem it is convenient to deal with Hilbert type structures, 
and therefore we assume in the further considerations that 17 £ i 2 (—f*, 0; "H). Since we do not 
assume the continuity of g in s £ [—f*, 0], we also need to add the (standard) initial conditions 
of the form u(t = 0) = uo(x) and d t u(t = 0) = -ui(x). 

The delayed system (without rotational inertia) is then given by: 

{ u tt + A 2 u + k 0 u t + f(u) + Lu = p 0 + g(u*,t) in flx(0,T), 
u = d u u = 0 on dClx(0,T), 
u( 0 ) = u o , u t (0) = u\, 
u |te(-t*,o) = V € L 2 (—t*, 0; Hq (fl)). 

Here the forcing term q(u l ,t) occurring on the RHS of the plate equation encompasses the gas 
flow. The operator L encompasses any spatial lower order terms which do not have gradient 
structure. 

We take a weak solution to EH on [0,T] to be a function 

u £ L 00 (0, T; H%(Sl)) n W£(0, T; L 2 (H)) n L 2 (-t*, 0; I/ 0 2 (fi)) 

such that the variational relation corresponding to (14.11) holds (see, e.g., [521 (4.1.39), p.211]). 
In order to consider the delayed system as a dynamical system with the phase space 

HeMx L 2 (n) x L 2 (-i*,0;«). 

With the above notation we introduce the operator S t : H e> H by the formula 


S t {u 0 ,u 1 , 77 ) = (u(t), «*(*),«*), 


(4.2) 


where u(t) solves (14.11) and u*(s) = u(t + s), s £ [—f*,0]. 

Theorem 4.1. Under the hypotheses imposed, EH generates a strongly continuous semigroup 
S t : H —> H which possesses a compact global and ’’smooth” attractor of finite fractal dimension O 

The main ingredient of the proof is a variant of hidden compactness of the term q(u*,t); 
indeed, a hidden regularity of the delay term leads to a variant of compensated compactness 
property which allows one to derive quasi stability inequality for the dynamics (see |24|1. We 
note that quasi-stability method originated in [20] (see also ED US [ 25 ], and is exposited in the 
recent monograph BSD- For delay models this method was also applied recently in [291 1 30] . 

Nonlinear PDEs with delay have been considered in various sources (see [72] and references 
therein). In relation to plate equations with delayed aerodynamic pressure, [22] (see also 151) 
provides a rather complete analysis of the delayed vK plate in the presence of rotational terms 
(and application to flow-plate interactions), and the references p3J [28] deal with the plate with 
delay in the presence of other nonlinear terms. One should stress that even at the formal 
level- with or without rotational inertia—the reduced delay problem is challenging. This is 
due to the fact that the underlying delay system is intrinsically non-gradient—the structure of 
the aeroelastic potential provides non-conservative and non-dissipative terms that destroy the 
Lyapunov nature of the energy relation. In view of this, results on asymptotic behavior are 
challenging. 

4 For the corresponding definitions and basic facts we refer to [221167 ]. 
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4.2 Piston Theory 

Classical piston theory (or law of plane sections 013108) replaces the acceleration potential of 
the flow, if = tr[(f> t + U<f> x ] with — [u t + Uu x \ in the RHS of the reduced plate equation (13.141) . 

Formally, we can arrive at the aforementioned, standard (linear) piston theory model by 
utilizing the reduction result in Theorem 13.61 a bound exists p.334]: 

l|g“(*)llo,n c H AM ( r )lln dr ( 4 - 3 ) 

for U sufficiently greater than 1. The constants c, C > 0 depend only on the diameter of f l. 
This indicates that the delay term, q u decays in the appropriate space as U increases. Hence, 
it is reasonable to guess that q u can be neglected in the case of large speeds U. @ So we arrive 
at the following model 


u tt + A 2 it + f(u) = po - [u t + Uu x ) in Slx(0,T), 
u = d v u = 0 on c>fix(0,T). 


This model was intensively studied in the literature for different types of boundary conditions 
and resistance damping forces, see [22] and also m for related second order abstract models. 
The typical result is the following assertion: 

Theorem 4.2. Under the conditions imposed on the model the equations m generates a 
dynamical system in the space i4g(f2) x possessing a compact global attractor of finite 

fractal dimension which is also “smooth”. 

A key avenue of research for supersonic flow-plate models in the standard panel configuration 
[2S| (see Section I3J1) is an explicit comparison between solutions to the piston theoretic plate 
model and solutions to the clamped flow-plate interaction as U becomes sufficiently large. This 
comparison should be made rigorous. This leads to the following fundamental question from 
modeling point of view: 

Open Question: Show that the associated solutions to (13.141) and (14.41) discussed above 
coincide for (i) t —> oo with U fixed and sufficiently large, or (ii) as U —> oo (the so called 
hypersonic limit) on some arbitrary time interval [Tj,T 2 ] (where T 2 is perhaps 00 ). Such 
results would provide rigorous justification (at the infinite dimensional level) for piston 
theoretic approaches to flow-plate models^ 

Remark 4.1. In addition to the classical, linear piston theoretic model discussed here (in (14.41) '). 
there are other piston theories. For instance, if we a priori assume that we deal with “low” 
frequencies regimes for the plate, then we can use the following expression for the aerodynamic 
pressure (see, e.g., mm and the references therein): 


p(x, t)=p Q - 


U 

Vu 2 - 1 


U 2 — 2 
U 2 - 1 


ut + Uu 



in the supersonic case (U >1). We note that this term has a fundamentally different structure, 
and its contribution to the dynamics is markedly different than the RHS of (14.41) . It is commonly 

5 See also Remark 1 1.1 1 for other (physically motivated) expressions the “piston” term). 

6 The reference m contains a result for the behavior of solutions to the piston theoretic plate (with RHS 
given by p* above) as U —> oo. The result, however, is only valid for arbitrarily small time intervals, and hence 
does not provide information about the behavior of solutions for arbitrary t. 
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recognized that the damping due to the aerodynamic flow is usually substantially larger in 
magnitude than the damping due to the plate structure. However as we can see the aerodynamic 
damping can be negative as well as positive. For instance, in the case 1 < U < \/2 we obtain 
negative damping |351 140| . This has a definite impact on the dynamics, producing instabilities 
in the system. 


5 Kutta-Joukowsky Boundary Conditions 

Arguably, the Kutta-Joukowsky boundary conditions for the flow (12.611 are the most important 
when modeling an airfoil immersed in a flow of gas mu. Not surprisingly, these boundary 
conditions —dynamic and mixed in nature—are also the most challenging from mathematical 
stand point. This was recently brought to evidence in a long treatise [2] . Various aspects of 
the problem—in both subsonic and supersonic regime—have been discussed in [2] within the 
context of mostly one dimensional structures. The aim of this section is to revisit the problem 
by putting it within the framework of modern harmonic analysis. The recently studies mm 
show how the flow condition (KJC) interacts with the clamped plate in subsonic flows in order 
to develop a suitable abstract theory for this particular model. Though the analysis is subsonic 
for (KJC), utilizing the flow energy from the supersonic panel (as above in Theorem 13.11) is 
effective in the abstract setup of this problem. In fact, even in the subsonic case, the analysis 
of semigroup generation proceeds through the technicalities developed earlier for supersonic 
case [26j . The key distinction from the analysis of the clamped flow-plate interaction (owing 
to the dynamic nature of the boundary conditions) is that a dynamic flow-to-Neumann map 
must be utilized to incorporate the boundary conditions into the abstract setup. The regularity 
properties of this map are critical in admitting techniques from abstract boundary control m , 
and are determined from the Zaremba elliptic problem |60] . To again decouple the dynamics 
into a dissipative piece and a perturbation, a trace regularity result must be established. The 
necessary trace regularity hinges upon the invertibility of an operator which is analogous (in two 
dimensions) to the finite Hilbert transform [2] |3) [68; . And this is a critical additional element of 
the challenging harmonic analysis brought about by the (KJC). When the problem is reduced in 
dimensionality to a beam structural model, this property can be demonstrated and our analysis 
has parallels with that in [21 [3] (and older references which are featured therein). Specifically, 
one must invert the finite Hilbert transform in L P (Q,) for 2; in higher dimensions, this brings 
about nontrivial (open) problems in harmonic analysis and the theory of singular integrals when 
(KJC) is in force. 

The full system with (12.61) can be written in terms of aeroelastic potential variable as: 


u tt + A 2 u + f(u ) = po + tr[ip] 


u = d v ii = 0 
(d t + Ud x )4> = ip 
{d t + Ud x )ip = A <p 
dv<p = -( d t + Ud x )u 
Sdt + Ud x )<j> = 0 


infix (0,T), 
on dfl x (0, T), 
in R+ x (0, T), 
in !hj_ x (0, T), 
onflx (0, T), 
on R 2 \H x (0, T), 


(5.1) 


with the given initial data (uq,Ui, (po,ipo). We make use of the flow acceleration multiplier 
(■ d t + Ud x )(f> = ip for the subsonic flow , taken with (KJC). Thus for the flow dynamics, instead 
of (</>•, cpt) we again have the state variables (</>•, ip) (as in the supersonic panel case) which then 
leads to a non dissipative energy balance, see (13.131) . 

Our result is formulated under the following regularity condition (to be discussed later). 
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Condition 1 (Flow Trace Regularity). Assume that 4>{x.,t) satisfies (12.41) and the Kutta- 
Joukowsky condition (KJC), and dttr[(j>\, d x tr[(j)\ £ £ 2 ( 0 , T; .ff _1 / 2_e (M 2 )) VT > 0. In 
addition, if = fit + Ufix satisfies the estimate 


IHV#)llff-i/2- e(R 2 ) dt < C T E f i{ 0) + / ||^(t)||n* 


(5.2) 


We note that the regularity required by Condition Q] is exactly the one that will make energy 
relation meaningful. To wit: u x (t) £ £Tq(I 2) and tr[if\ £ L 2 (0, T; H~ e ), 9 > 1/2 defines the 
correct duality pairing in the tangential direction on M 2 . It is also at this point where we use 
the fact that u satisfies it = 0 on dfl . Thus simply supported and clamped boundary conditions 
imposed on the structure fully cooperate with this regularity. The principal result of [52] (see 
also E2) reads as follows: 

Theorem 5.1. With reference to the model & with 0 < U <1: Assuming the trace regularity 
Condition^ holds for the aeroelastic potential if = {(f t + U(f x ), there exists a unique finite energy 
solution which exists on any [0,T], This is to say, for any T > 0, (<j>, (f> t , it, u t ) £ C(0,T;Y) for 
all initial data (</> 0 , </>i, uq, iq) £ Y. This solution depends continuously on the initial data. 


The proof of Theorem 15.11 follows the technology developed for the supersonic case in [26] 
and is given in |52j . In view of Theorem 15.11 the final conclusion on generation of nonlinear 
semigroup is pending upon verification of the flow regularity Condition. While the complete 
solution to this question is still unavailable, and pending further progress in the theory of finite 
Riesz transforms, we can provide positive answer in the case when 12 is one dimensional. In fact, 
this positive assertion also follows a-posteriori from direct calculations given in PQ, and based 
on the analysis of finite Hilbert transforms. The arguments given below are independent and 
more general, with potential adaptability to multidimensional cases. The corresponding result 
is formulated below. 

Theorem 5.2. Assume 12 = (—1,1) (suppressing the span variable y). Then the flow trace 
regularity in Condition [7] holds for <j>. In this case the semiflow defined by 15.21 ) taken with 
(KJC) generates a continuous semigroup. 

In what follows below we shall address validity of Condition [T] in the one dimensional case. 
In particular, we wish to provide the relation to L p , p / 2, invertibility of Hilbert/Riesz 
transforms^ The thrust of the argument is the trace regularity of the following flow problem in 


M-* 


' (d t + Ud x ) 2 (j) = Afi 

</(0) = </ 0 ; </t(0) = cf>i, 

d v <f = /i(x, t ) 

tr[(ft + U<j) x \ = tr[ip\ =0, x £ R 2 \I2 


in R^_ x (0,T), 
on 12 x (0, T) 


(5.3) 


where 2i(x, t), defined on 12, is the downwash generated on the structure and tr[if\ is the aeroe¬ 
lastic potential. The flow trace regularity problem can be rephrased as the following inverse 
type boundary-boundary problem: 

Given h(x,y,t) £ L 2 (K x 12), find if, and hence tr[ip\, on 12. 

In line with the similar analyses in |561159] , we take zero initial data (later we may apply 
the principle of superposition). We consider the Fourier-Laplace transform of the original linear 


7 More complete analysis and details are given in EH. 
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equation (formally, sending (x,y) —> i{j] x ,rj y ) and t —> r = (a + i/3): 


(r + iUr] x ) 2 (j) + \r)\ 2 (f) = (t 2 + 2Uir) x r + (1 - U 2 )r] 2 x + rfyj) = d 2 (p. 
Let us introduce the following PDO symbol: 

D(v, t) = (t 2 + 2 Uif] x r + (1 - U 2 )r] 2 + ? f y ), 
and so we must solve (in z) the following ODE: 

d z 2 (j> = D(T,rj)4>, 


whose general solution is 


< K t ’V,z ) = <j>(r,r], 0)e 


— zy/D 


Our next step is to relate the boundary conditions to the the normal derivative of and the 
acceleration potential ip: 

V>(D V, z) = (t + iUr) x )<j>(T, rj, 0)e~ zy/B 
d z cp = 0) 


(5.4) 


Hence, on we have the relation 


d z (p = —ip = m(r, rj)ip 


with 


t + iUr] x 


l (Vi T ) = 


-Vd 


(5.5) 


t + iUr] x 

which corresponds to zero order PDO £ S' 0,1 j43( [66]. This relation is supplemented with the 
information 


d z (p\ z =o = h € L 2 (R x n) (5.6) 

ip(t,x,y, 0) = 0, outside D. (5.7) 

Let & denotes the operator corresponding to the multiplier m(r],T) in (15.41) . Let Pq : R 2 — > Q 
denote projection on Q. Then for any ip which is zero off D, the problem stated above can be 
recast as: given h , find ip so that 


h = Pn^ip, x £ D, t > 0 (5-8) 

noting that & acts on a “truncated” function ip, so that the PDO Pq3? can be viewed as PDO 
on x (0, oo) (viewing L 2 {^L) as a subset of L 2 (R 2 ) with associated projection). We construct 
3" as follows: for any ip with ip = 0 outside D, we construct Laplace (time) Fourier (space) 
transform of ip which we denote ip{T,rj). Then the operator 3 is the PDO such that 

{3ip)(T,r]) = m(T,T])ip(T,T]) (5.9) 

Taking the inverse of the Fourier-Laplace transform, and applying projection on Q gives the 
actual definition of the operator. Thus, the problem of trace regularity (PDO Neumann-Dirichlet 
map) reduces to solving (15.81) with the PDO given by (|5.9I) and Mikhlin type of multiplier m(r, rj). 
The relationship between ip and h depends on the properties of the Mikhlin multiplier to ( t , rf) 
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in the class of symbols S' 0,1 43], [66]. In fact, as we will show below, & is related to the inverse of 
a finite Hilbert transform in the x variable. We make this notion precise below, but we do note 
at this stage that the integral transform which arises is not standard in the two dimensional 
scenario. It is only when when Q is reduced to an interval that the multiplier above contains the 
classical finite Hilbert transform EH EH]. Otherwise we deal with a “finite” Riesz transform. If 
we denote by H. Hilbert transform with the symbol symb{T~L ) = —i-^-and by Eq the extension 
operator (by zero) from L p {fX) onto L p (R”), then finite Hilbert transform can be written as 

H f = PnUE n e ->• L p (fi)), 1 < p < oo. 


When Q = I this corresponds to a classical finite Hilbert transform. For higher dimensions one 
can consider Q as an union of intervals -generalizing the concept to finite Riesz transforms. The 
main point critical to the analysis is the fact that m 

Hf 1 € Sf(L p (Sl)),p < 2 (5.10) 

Remark 5.1. The connection between integral equations appearing in the study of aeroelasticity 
and invertibility of finite Hilbert transforms has been known for many years HH, and dates back 
to Tricomi. This approach has been critically used in 01 where the analysis is centered on 
solvability of integral equations of Possio’s type connecting the downwash with the aeroelastic 
potential. We follow the same conceptual route as these references, with different technical tools. 
Our approach is based on microlocal analysis—rather than explicit solvers of integral equations 
arising in a very special case (of a purely one dimensional setting). Though our final estimate 
depends on an assumption which has only been demonstrated for the one dimensional case, 
we believe that the microlocal approach provides new ground for extending the flow-structure 
analysis to multidimensional settings. 

In line with hyperbolic character of the symbol, it is convenient to introduce the variable 
Zjj = + U. In the case when the dimension of H is equal to one our multiplier then becomes, 


m(z,r]) 


V 


r(z,rj) 


m 0 (v)r(z,v) 


where 


with 


r{z,rj) 


VD~i 

Zjj — iar)- 1 


D\{z, rf) = —a 2 ri 2 - 1 + z% - 2i — (zu). 

V 


We are interested in the invertibility (and regularity properties of the inverse) of the operator 
associated to the multiplier m{z,rf). Noting that mo(r]) is a symbol corresponding to Hilbert 
transform, we are seeking a bound from below for r[z,rf). The critical regions for the analysis 
of this symbol occur for large values: \r]\ >> R , |/?| >> R , for R sufficiently large. Analyzing 
the symbol r(z,r]), we note that when \z\ —> 0 and |^| —> oo, the symbol r(z,rj) converges to 
constant, nonzero values. A loss of regularity occurs when |^| —>• 1. In that case, the asymptotic 
behavior is characterized by 



This indicates a loss of 1/2 of a spatial derivative, i.e. a bound from below of the form 


r(z, r])y/rj > c > 0. 

The above leads to the following lemma: 
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Lemma 5.3. For Q = I C R a compact interval, the operator associated with the multiplier 
m{rf) is invertible from L p (fl) —> W~ 1 ^ 2 ,p {yi) for p £ (1,2). 

Proof. We first reference the fact that finite Hilbert transform is Fredholm on L p {I) for p £ 
(1,2) U (2,oo), in particular (15.101) takes place. 22- Moreover, the operator associated with 
m o{ T l) T l ~ 1 maps W _1 ’ p (/) onto L p {I ), for the same range p £ (1,2) U (2, oo). This follows from 
derivative formula for finite Hilbert transform for functions vanishing outside I [71] . Thus, by 
interpolation, the operator associated to the multiplier m(z,r/) is invertible from L P (I) into 
W~ 1 ' 2,P (I). Since our analysis is uniform in ft, we have that for h £ L2((0, T ) x /)), one obtains 
that aeroelastic potential tr[if>\ is in L2((0, T), W _1 / 2,p (/)) for p < 2. □ 

By Sobolev’s embeddings W _1 / 2 ’ p (/)) C with e(p) > 0 possibly arbitrarily 

small -justifying the regularity of tr[ip\ in Condition [T| when ft = I. 

Remark 5.2. As discussed above, the generation of semigroups for an arbitrary three dimensional 
flow is subjected to the validity of the trace Condition [T] While it is believed that this property 
should be generically true, at the present stage this appears to be an open question in the 
analysis of singular integrals and depends critically on the geometry of f l in two dimensions. 

Remark 5.3. We note that the regularity of aeroelastic potential required by Condition Q] in one 
dimension, follows from the analysis in Eiia, where the author proves that aeroelastic potential 
ip £ £ 2 ( 0 , T, L q (Q)) for q < 4/3 . Since for p > 4 there exists e > 0 such that 

H 1 / 2+e (H) C L p (Q), p > 4, dim Q < 2, 

and one then obtains that L q {kV) C 77 _1 / 2_e (H) with q < 4/3. 

Remark 5.4. The loss of 1/2 derivative in the characteristic region was already observed and 
used in the analysis of regularity of the aeroelastic potential for the Neumann problem with 
supersonic velocities U ES- However, in the case of (KJC) there is an additional loss, due to 
the necessity of inverting finite Hilbert transform which forces to work with L p theory for p < 2. 
This is due to the fact that finite Hilbert transform is invertible on L p , p < 2 -rather than p = 2. 

The author of mm considers a linear wing immersed in a subsonic flow; the wing is taken 
to have a high aspect ratio thereby allowing for the suppression of the span variable, and 
reducing the analysis to individual chords normal to the span. By reducing the problem to a 
one dimensional analysis, many technical hangups are avoided, and Fourier-Laplace analysis is 
greatly simplified. 


6 Some Experimental and Computational Observations 
Reinforced by Mathematical Analysis 

In what follows we shall present several observations which, on one hand confirm experimentally 
and numerically some of the mathematical findings discussed before and, on the other hand, 
raise open questions and indicate new avenues for mathematical research. 

1. Stability induced by the flow in reducing dynamics to finite dimensions: Exper¬ 
imentally we see that the flow has the ability of inducing stability in the moving structure. 
This is the case when the structure itself does not possess any mechanical or frictional 
damping. This dissipative effect is not immediately noticeable in the standard energy 
balance equation. However, rigorous analysis reveals that the stabilizing effect of the flow 
reduces the structural dynamics to a finite dimensional setting. It is remarkable that 
such conclusion is obtained from mathematical considerations—as it is not clearly arrived 
at by either numerics or experiment. 
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2. Non-uniqueness of final nonlinear state: Our analysis of asymptotic behavior asymp¬ 
totically reduces the panel dynamics (in the subsonic case) to the set of stationary solu¬ 
tions. However, such stationary solutions are not unique in general |f| In this case it 
is possible that there exist several locally stable equilibria in the global attractor. This 
makes it possible to explain the fact—observed experimentally—that the buckled plate in 
an aerodynamic flow does not have a final, unique nonlinear state; it can depend on the 
path which is traveled in the parameter space, and/or initial conditions. See [401 and 
references cited therein, as well as [HEMS]. The latter discuss the chaotic responses 
that may occur for a buckled plate in a gas flow. 

3. Surprisingly subtle effects of boundary conditions: As already noted, as the struc¬ 
ture’s boundary conditions are changed, so is the dynamic stability/instability of the 
system. Experimental and numerical studies reveal that a change of boundary conditions 
(from clamped to free plates, for instance) can lead from flutter (resulting in periodic 
solutions), to divergence (resulting in a static buckling). Computationally and experimen¬ 
tally it is found that in subsonic flows the trailing edge boundary condition is the key. If 
the trailing edge is fixed, then divergence will occur. If it is free, then flutter will occur. 
A recent computational and experimental study of the effects of boundary conditions on 
plate flutter, divergence, and limit cycle oscillations has been described in [47] . Our the¬ 
oretical results on convergence to equilibria confirm these experimental studies—at least 
with respect to fixed (clamped) trailing edge. 

4. Limitations of VK theory and new nonlinear plate theory: When the leading edge 
is clamped and the side edges and trailing edge are free flutter and periodic oscillations 
have amplitudes on the order of the plate chord. Thus vK theory is no longer accurate. 
This may explain our mathematical difficulties in carrying the vK analysis in configurations 
where a portion of the boundary is free. However, a novel, improved nonlinear plate theory 
has been developed and explored computationally, and correlated with experiment. This 
provides a novel, challenging, and exciting opportunity for mathematical analysis. A new 
nonlinear structural theory for cantilevered plates has been developed in [65] and applied 
to aeroelastic studies including a comparison with experiment in [63] . 

5. Aerodynamic theory: Very recently piston theory has been extended to lower supersonic 
Mach numbers. This may be helpful to prove mathematically some of the results in 
the spirit of those given above [Ml HO]- Indeed, our analysis of delayed plate indicates 
that, ultimately, the delayed term (retained for subsonic velocities) does not have any 
major effect on qualitative behavior. It is rather the dispersive effect (emanating from the 
— Uu x term that leads to potential instabilities in the dynamics. 

6. Extensions to incorporate solar radiation: The pressure due to solar radiation has a 
similar mathematical form to that of piston theory and has been of recent interest in the 
context of interplanetary transportation using solar sails. Such sails may also experience 
flutter due to solar pressure. The nonlinear models for pressure differ, however (aero¬ 
dynamic versus solar radiation), and represent an opportunity for further mathematical 
studies. See the work of Dowell in [M] and Gibbs and Dowell in [45] . 

8 This can be seen from the results in Section 6.3.6 in m in the subsonic case. 
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7 Open Problems 

In what follows below we shall list some open problems which naturally emerge from the material 
presented above. 

1. Configurations involving a free-clamped plate: As pointed out in section 1.3, the 
free-clamped boundary conditions imposed on the plate are of great physical interest. 
From the mathematical point of view, the difficulty arises at the level of linear theory 
when one attempts to construct ” smooth” solutions of the corresponding evolution. Typ¬ 
ical procedure of extending plate solutions by zero outside Q leads to the jump in the 
Neumann boundary conditions imposed on the flow. In order to contend with this is¬ 
sue, regularization procedures are needed in order obtain smooth approximations of the 
original solutions. While some regularizations have been already introduced in [26] , more 
studies are needed in order to demonstrate the effectiveness of this ’’smoothing” for the 
large array of problems described in this work. Free boundary conditions in the context 
of piston theory and boundary dissipation have been recently studied in [8] and also [22] • 

2. Kutta-Joukowsky boundary conditions for the flow in higher dimensions: As 

explained above, the (KJC) are used in the models of panels that are partially free. Jump 
of the pressure off the wing is a typical configuration [2]. From the mathematical point of 
view the difficulty lies, again, at the level of the linear theory. In order to deal with the 
effects of the unbounded traces tr[ip\ in the energy relation microlocal calculus is necessary. 
This has been successfully accomplished in [25] where clamped boundary conditions in 
the supersonic case were considered. However, in the case of (KJC) there is an additional 
difficulty that involves “invertibility” of finite Hilbert (resp. Riesz) transforms. This latter 
property is known to fail within L 2 framework, thus it is necessary to build the L p theory, 
p / 2. This was for the first time observed in [3] and successfully resolved in the one 
dimensional case. However, any progress to higher dimensions depends on the validity of 
the corresponding harmonic analysis result developed for finite Riesz transforms. Some 
preliminary results in the subsonic case are in [52] . 

3. The transonic regime: The treatment presented above excludes the transonic velocity 
U = 1. Indeed, for U = 1 the analysis provided breaks down in the essential way [7]. 
While numerical, experimental work predicts appearance of shock waves [36] . to our best 
knowledge no mathematical treatment of this problem is available at present. 

4. Nonlinear flow models: Finally, the ultimate goal is to consider a fully nonlinear flow 
model. Experimental-numerical results predicting shock waves in the evolution are par¬ 
tially available [35]. However the mathematical aspects of this problem are presently wide 
open. 

5. Rigorous comparison (asymptotics) of delay plate dynamics and piston the¬ 
oretic dynamics: There are also a question of justification of the “piston” theory at 
the level of global attractors via hypersonic limit U —> + 00 . This requires asymptotic 
(U —> 00 ) closeness of the corresponding solution. (See the Open Question earlier in the 
text.) 
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